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§ Abstract 

The COMPASS experiment at the CERN SPS has studied the diffractive dissoci- 
ation of negative pions into the ir~ir~ir + final state using a 190GeV/c pion beam 
hitting a lead target. A partial wave analysis has been performed on a sample 
of 420 000 events taken at values of the squared 4-momentum transfer t' between 
0.1 and lGeV 2 /c 2 . The well-known resonances ai(1260), a2(1320), and 7^(1670) 
are clearly observed. In addition, the data show a significant natural parity ex- 
change production of a resonance with spin-exotic quantum numbers J PC = 1 
at 1.66GeV/c 2 decaying to pn. The resonant nature of this wave is evident from 
the mass-dependent phase differences to the J PC = 2 h and 1 ++ waves. From 
a mass-dependent fit a resonance mass of (1660 ± lOigj MeV/c 2 and a width of 
(269 ± 21±f 4 ) MeV/c 2 are deduced, with an intensity of (1.7 ± 0.2)% of the total 
^ intensity. 
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In the SU(3)fl avor constituent quark model, light mesons are described as bound 
states of a quark q and an antiquark <f with quark flavors u, d, s. Mesons are classified 
in J PC multiplets, with the total angular momentum J, the parity P, and the particle- 
antiparticle conjugation parity C, which is defined through the neutral flavorless members 
of a given multiplet. The isospin / and the G-parity further characterize mesons containing 
light quarks. In the quark model, P, C and G are given by 

p=(_l) L + 1 ) C=(-1) L+S , G=(-1) I+L+S , (1) 

where L is the relative orbital angular momentum of q and q', and S the total intrinsic 
spin of the qq' pair, with S = 0,1. The constituent quark model has been quite successful 
in explaining many of the properties of mesons as well as, to a large extent, the observed 
meson spectrum, even though it makes no assumptions concerning the nature of the 
binding force, except that hadrons are postulated to be color-singlet states. In Quantum 
Chromodynamics (QCD), the interaction between colored quarks is described by the 
exchange of gluons which carry color themselves. Owing to this particular structure of 
QCD, color-singlet mesons can be formed not only by constituent quarks, but also by other 
configurations like four-quark objects or gluonic excitations. These non-gg' configurations, 
however, will mix with ordinary qq' states with the same J PC , making it difficult to 
disentangle the contribution of each configuration. The observation of exotic states with 
quantum numbers not allowed in the simple quark model, e.g. J PC = , + ~, 1 h , . . ., 
would give clear evidence that quark-gluon configurations beyond the quark model, as 
allowed by QCD, are realized in nature. 

The lowest-lying hybrid, i. e. a system consisting of a color octet qq' pair neutralized 
in color by a gluonic excitation, is expected [1] to have exotic quantum numbers J PC = 
1 h , and thus will not mix with ordinary mesons. Its mass is predicted in the region 
1.3 — 2.2 GeV/c 2 . The systematics of hybrid meson production and decay has been worked 
out in the flux-tube model [2]. There are three experimental candidates for a light 1~ + 
hybrid. The 7Ti(1400) was observed by E852 [3] and by VES [I] in the reaction n~N — > 
i]n~N, and by Crystal Barrel [5j [6] in pn — > n~n°r] and pp — > 2n°r] Dalitz plot analyses. 
Another 1~+ state, the tti(1600), decaying into pn [TJ El E] , rfir PH [□] , /i(1285)tt [HH3], 
and &i(1235)7r [T31 H3] was observed in peripheral n~p interactions in E852 and VES, and 
confirmed in pp — > binn [15]. The resonant nature of both states, however, is still heavily 
disputed in the community [U [13] . In a different analysis of a larger data set of E852 no 
evidence for an exotic resonance at 1.6 GeV/c 2 in the 3n final state was found [16J. A third 
exotic state, 7^(2000), decaying to fin and bin, was seen in only one experiment [T2"j [T4]. 

In order to shed new light on these questions, the COMPASS collaboration, operat- 
ing a large-acceptance and high-resolution spectrometer [TF] situated at the CERN Super 
Proton Synchrotron (SPS), is gathering high-statistics event samples of diffractive reac- 
tions of hadronic probes into final states containing both charged and neutral particles. 
Diffractive dissociation is a reaction of the type a + b —> c + d with c— > 1 + 2 + ■ ■ ■ + n, 
where a is the incoming beam particle, b the target, c the diffractively produced object 
decaying into n particles, and d the target recoil particle, with 4-momenta p a . . . pd, re- 
spectively. The production kinematics is described by two variables: s and t' = \t\ — |t| m in, 
where s = (p a + pb) 2 is the square of the total center of mass energy, t = (p a — p c ) 2 is 
the square of the four momentum transferred from the incoming beam to the outgoing 
system c, and |t| m i n is the minimum value of \t\ which is allowed by kinematics for a given 
mass m c . 
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Figure 1: Invariant mass of the 3n system for 0.1 GeV 2 /c 2 < t' < 1.0 GeV 2 /c 2 (histogram), 
and intensity of the background wave with a flat distribution in 3-body phase space 
(triangles), obtained from a partial wave analysis in 40MeV/c 2 bins of the 3tt mass and 
rescaled to the binning of the histogram. Both the invariant mass spectrum and the 
background distribution are not acceptance corrected. 

First studies of diffractive reactions of 190GeV/c tt~~ on a 3 mm lead target were 
carried out by COMPASS in 2004. The 7r~7r~7r + final state was chosen because the dis- 
puted 71"! (1600) meson with exotic J PC had previously been reported in this channel. The 
trigger selected events with one incoming particle and at least two outgoing charged par- 
ticles. In the offline analysis, a primary vertex inside the target with 3 outgoing charged 
particles is required. Since the recoil particle was not detected, the following procedure is 
applied in order to select exclusive events. The beam energy E a is very well approximated 
by the measured total energy E c of the 3n system with a small correction arising from the 
target recoil, which can be calculated from the measured scattering angle 9 = l{p a ,Pc), 
assuming that the target particle remained intact throughout the scattering process. Then 
an exclusivity cut is applied, requiring E a to be within ±4 GeV of the mean beam energy. 
Events with a wide range of t' from zero up to a few GeV 2 /c 2 were recorded. For the anal- 
ysis presented in this letter we restrict ourselves to the range where candidates for spin 
exotic states have been reported in the past: 0.1 GeV 2 /c 2 < t' < 1.0 GeV 2 /c 2 , far beyond 
the region of coherent scattering on the Pb nucleus. Figure [T] shows the invariant mass of 
the corresponding events. In our sample of 420 000 events in the mass range between 0.5 
and 2.5GeV/c 2 , the well-known resonances ai(1260), a 2 (1320), and 7r 2 (1670) are clearly 
visible in the 3ir mass spectrum. 

A partial wave analysis (PWA) of this data set was performed using a program 
which was originally developed at Illinois [IB], and later modified at Protvino and Mu- 
nich. An independent cross-check of the results was performed using a different PWA 
program developed at Brookhaven [19] and adapted for COMPASS [20]. At high a/s, the 
reaction can be assumed to proceed via t-channel Reggeon exchange, thus justifying the 
factorization of the total cross section into a resonance and a recoil vertex without final 
state interaction. The exchanged Reggeon may excite the incident pion (J p = 0~) to a 
state X with different J p , limited only by conservation laws for strong interactions. For 
the (37r) - final state / > 1; we assume 1 = 1 since no flavor-exotic mesons have been 
found. Since in addition G = — 1 for a system with an odd number of pions, C = +1 
follows from eq. [T] We take the phenomenological approach of the isobar model, in which 
all multi-particle final states can be described by sequential two-body decays into inter- 
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mediate resonances (isobars), which eventually decay into the final state observed in the 
experiment. All known isovector and isoscalar tctc resonances have been included in our 
fit: (tctt)s (comprising the broad a(600) and /„(1370)), p(770), /o(980), / 2 (1270), and 
p 3 (1690) [Sj. It is possible that there exists a direct 3-body decay into (3tc)~ without an 
intermediate di-pion resonance; in the isobar model, such a decay mode without angular 
correlations is represented by <r(600)+7r~ with L = and J p = 0" . Possible complications 
to the isobar model from unitarity constraints are not an issue here; such effects enter in 
the formulation of the model only when all possible decay modes are simultaneously fit, 
which may include the final states containing 7r°, rj, rf , u , KK or NN . The spin-parity 
composition of the excited state X is studied in the Gottfried- Jackson frame, which is 
the center of mass frame of X with the z-axis along the beam direction, and the y-axis 
perpendicular to the production plane, formed by the momentum vectors of the target 
and the recoil particle. 

The P WA is done in two steps. In the first step, a fit of the probability density in 
37r phase space is performed in 40MeV/c 2 bins of the 3tt invariant mass m (fit in mass 
bins). No dependence of the production strength for a given wave on the mass of the 3n 
system is introduced at this point: 



N r 



a indep (T,m,t') = 

e=±l r=l 



m: 



N{(m) = J \^(r',m)\ 2 dr' 



(2) 

Here, Tf r are the production amplitudes and ipl the decay amplitudes, the indices % and e 
denoting different partial waves, characterized by a set of quantum numbers J PC M e [isobar] L; 
M is the absolute value of the spin projection onto the z-axis; e is the reflectivity [21], 
which describes the symmetry under a reflection through the production plane, and which 
is defined such that it corresponds to the naturality of the exchanged Regge trajectory; L 
is the orbital angular momentum between the isobar and the bachelor pion. The different 
t' dependence of the cross section for M = and M = 1 states is taken into account by 
including different functions of t', flit') oc exp (— bt') (M = 0) and flit') oc t'exp(— bt') 
(M = 1), where the slope b has been obtained from the data by first making fits in slices of 
t'. The 3-body decay amplitudes are constructed using non-relativistic Zemach tensors 
[22] . They are properly Bose-symmetrized to take into account the combinatorics due to 
the spin-0 nature of the final state pions. They depend on the set of five parameters r 
specifying the 3-body decay kinematics, but do not contain any free parameters. The nor- 
malization factors Nf(m) contain angular- momentum barrier factors and quasi- 2-body 
phase space factors, taking into account the non-zero widths of isobars. Dividing each 
decay amplitude by N^(m) compensates its dependence on the mass inside each mass 
bin. Equation [2] includes a coherent sum over waves with different J PC M, allowing them 
to interfere. It also contains two non-coherent sums over the reflectivity e and the rank 
N r [21]. Assuming that the recoiling target particle is a nucleon, and neglecting nuclear 
effects, we set N r = 2, corresponding to a mixture of helicity-flip and helicity-non-flip 
processes at the baryon vertex. A total of 42 partial waves are included in the first step of 
the fit. It comprises the non-exotic positive-reflectivity waves with J PC = h (M = 0), 
l ++ ,2-+,3 ++ ,4- + (M = 0,1), 2 ++ ,4 ++ (M = 1), the exotic 1~+ (M = 1), and the 
negative-reflectivity waves 1 h ,2 ++ (M = 0,1), l ++ ,2 h (M = 1), taking into account 
all relevant known decay modes into the isobars listed above. The fit also contains a 
background wave, characterized by a uniform distribution in 3-body phase space, which 
is added incoherently to the other waves. In the fit the elements of the spin-density ma- 
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trix, given by /?f- = J2 r Ti r Tj*, are determined simultaneously for all 42 waves in a given 
mass bin using an extended maximum likelihood method. The diagonal elements of the 
spin-density matrix are the intensities of the corresponding waves, while the off-diagonal 
elements determine the phase differences between two waves. There is no a-priori assump- 
tion on a resonant behavior of a given wave at this first step. The fit also takes into account 
the experimental acceptance of the spectrometer, calculated from a phase-space Monte 
Carlo simulation of the apparatus. It is worth stressing that COMPASS has an excellently 
uniform acceptance for diffractively produced 3n events of the order of 60% over the whole 
phase space. In order to verify that indeed the global maximum has been found by the 
fit, up to 100 attempts with randomly chosen start parameters are performed for each 
mass bin. If multiple solutions are found within one unit of log likelihood, the average of 
the two extreme solutions for each parameter is used. The difference of the two extreme 
solutions is added linearly to the statistical error of the best fit for each parameter. 

In the second step of the PWA a \ 2 fit of the spin-density matrix elements obtained 
for each mass bin in the first step is performed in the mass range from 0.8 to 2.32 GeV/c 2 , 
taking into account the mass dependence of the produced resonances (mass-dependent fit). 
The elements of the spin density matrix are expressed as p\^ = J2 r Al r (m)Aj*(m) with 
amplitudes A\ r {m) = J2kCi kr BWk(m) Nf(m), with BWk(m) denoting relativistic Breit- 
Wigner functions (with constant or dynamic widths depending on whether branching 
ratios of the corresponding resonance are known) or, where required by the fit, a coherent 
background reflecting non-resonant production of the corresponding partial wave, e.g. via 
the Deck effect [23J. For the latter, an empirical parameterization consisting of a simple 
exponential, exp(— ap 2 ), with p being the break-up momentum for the 2-body decay of the 
produced resonance and a a fit parameter is used. In the mass-dependent fit the complex 
production amplitudes Cf kr) and the parameters of BWk(m) are determined for a subset 
of six waves, the selected waves showing either significant amplitudes or rapid relative 
phase changes in the 1.7 GeV/c 2 mass range: 0" + + / (980)vr S, 1++0+ pix S, 2~ + + / 2 tt S, 
2++1+ pnD, 4 ++ 1 + pttG, and the exotic l^+l+pyrP. 

The intensity of the background wave resulting from the fit in mass bins is included 
in Figure [l] (triangles). Figures |2]^a)-(c) show the acceptance-corrected intensities of the 
three most prominent waves l ++ + pnS, 2~ + + f 2 n S, and 2 ++ l + pn D, respectively, as 
determined from the fit in mass bins (data points with error bars). The data also show 
a significant natural parity exchange production of a wave with spin-exotic quantum 
numbers J PC = 1 h at 1.66GeV/c 2 decaying to pn (P wave), presented in Fig. [2^d). The 
resonant nature of the exotic wave is evident from its phase differences to the dominant 
l ++ + pn S and 2~ + + f 2 ^ S waves, shown in Figs.[3|a) and (b), respectively (data points). 
For the latter, shown in Fig.J^b), no significant change in the phase difference between 1.4 
and 1.9 GeV/c 2 is observed, which is attributed to the fact that there are two resonances, 
7Ti(1600) and 7^(1670), with very similar masses and widths, causing the relative phase 
difference to be almost constant. In contrast to this the phase difference to the 1 ++ 
wave, shown in Fig. [3^a), clearly shows an increase around 1.7GeV/c 2 . As the ai(1260) 
is no longer resonating at this mass, this observation can be regarded as an independent 
verification of the resonating nature of the 1 _+ wave. 

The solid lines in Fig. [2] show the total intensity from the mass-dependent fit for the 
corresponding waves. For the l ++ + pn S wave shown in Fig. |2](a) it is well known that 
there is a significant contribution of non- resonant production through the Deck effect [21] , 
indicated by the dotted line. Its interference with the ai(1260) (dashed line) shifts the 
peak in the data to a slightly lower value than the peak position of the resonance. The 



4 




Figure 2: Intensities of major waves l ++ + pn S (a), 2~ + + /2 71 " S (b), and 2 ++ l + pn D (c), 
as well as the intensity of the exotic wave 1 h l + pit P (d), as determined in the fit in mass 
bins (data points with error bars). The lines represent the result of the mass-dependent 
fit (see text). 



2 h + /27r S wave shown in Fig. |2[b) is well described by a single resonance, the 7^(1670). 
The 2 ++ l + pit D wave displayed in Fig. [2|c) is dominated by the a2(1320) with a small 
contribution from the 02 (1700), whose parameters have been fixed to PDG values be- 
cause of the limited statistics. The intensity of the exotic 1~ + 1 + pn P, shown in Fig. [2^d), 
is well described by a Breit-Wigner resonance with constant width at 1.66 GeV/c 2 (dashed 
line), which we interpret as the 7Ti(1600), and a non-resonant background (dotted line) at 
lower masses. The resonant component of the exotic wave is strongly constrained by the 
mass-dependent phase differences to the l ++ + pn S and the 2 h + f^R S waves, which are 
well reproduced in the mass-dependent fit (solid lines in Fig. [3]). The parameters deduced 
for the masses, widths and intensities of the resonances included in the mass-dependent 
fit are given in Table [TJ where the first uncertainty corresponds to the statistical error, 
the second to the systematic error. The intensities are given for the resonant part of the 
corresponding wave integrated over the mass range from 0.8 to 2.32 GeV/c 2 , and are nor- 
malized to the total intensity from the mass-dependent fit, corresponding to 38.7(2)% 
of the acceptance-corrected data sample in the same mass range. The intensity of the 
7Ti(1600) is found to be (1.7 ± 0.2)% of the total intensity, while the non-resonant contri- 
bution is (1.5 ± 0.2)%, with a linear correlation coefficient between the two intensities of 
15%. The dominance of natural- over unnatural-parity exchange is more pronounced than 
in the BNL case at 18GeV/c [S]. This is possibly due to the decreasing contribution of 
unnatural-parity exchange with an increasing beam energy, if the natural-parity exchange 
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Figure 3: Phase differences of the exotic 1 h l + pn P wave to the l ++ + pre S (a) and the 
2~ + + /2 71 " S (b) waves. The data points represent the result of the fit in mass bins, the 
lines are the result of the mass-dependent fit. 
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Table 1: Resonance masses, total widths, and intensities for the specified decay channel 
of the six waves included in the mass-dependent fit to the data. The first uncertainty 
corresponds to the statistical error, the asymmetric second one to the systematic error. 
The last two columns give the corresponding PDG values [26] . 



is mediated mostly by the Pomeron. 

The systematic errors were estimated from the data by testing the stability of the 
result with respect to various assumptions made in the analysis, e.g. adding or removing 
certain waves, varying cuts or initial parameters for the fit. One such study concerns the 
choice of the rank N r used in the PWA. Although N r = 2 is physically motivated from 
the fact that, at high t', incoherent diffraction from individual nucleons dominates the 
reaction, fits with N r = 1 and 3 were tried as well. The intensity in the background 
wave relative to the total acceptance-corrected data sample in the mass range from 0.5 
to 2.5GeV/c 2 increases from 5.8% for N r = 2 to 19% for A^ r = 1, while it drops to 
1.2% for iV r = 3. At the same time, however, N r = 3 was found to cause larger bin-to- 
bin fluctuations without significantly altering the result. Given the level of the present 
statistics, we therefore conclude that the optimum rank is N r = 2. In an attempt to 
account for the low-mass shoulder in the intensity of the l~ + l + p7rP wave we also tried 
to include a 7Ti(1400) into the mass-dependent fit, with parameters fixed to PDG values 
[25] . This reduced the background intensity to a negligible value and shifted the resonance 
mass of the 7Ti(1600) to a slightly smaller value, which is reflected in its systematic error, 
but did not affect the intensity or the phase differences of any of the other waves in 
the mass-dependent fit. Releasing the parameters of the 7Ti(1400), however, causes the 
fit to become unstable. This can be attributed to the fact that the 7Ti(1400), if present 
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at all, couples only weakly to the 7r~7r~7r + final state. Another study of the sensitivity 
of the 7Ti(1600) intensity to the functional form of the background was performed using 
a background parameterization without angular-momentum barrier factors, which did 
not alter the result. Other systematic studies included a shift of the 40MeV/c 2 mass 
bins by 20MeV/c 2 , the use of rotation functions with relativistic factors [2Z] instead of 
Zemach tensors for the fit in mass bins, and the inclusion of four waves with M — 2. 
The use of different parameterizations for the a and p mesons also did not influence the 
result. Performing the fit in mass bins on non-exclusive events, i.e. events outside the 
aforementioned exclusivity cut, no signal is observed in the 1~ + wave. 

An incomplete acceptance of the spectrometer, not properly taken into account in 
the Monte Carlo simulation, or an incomplete set of waves may introduce leakage of non- 
exotic waves into the 1 h wave. In order to study this effect, Monte Carlo events were 
generated using the parameters of 16 dominant waves, excluding the which were de- 
termined in a mass-dependent fit, and simulating the decay patterns of the corresponding 
decay channels. Performing the same PWA for the Monte Carlo data as for the real data 
it was found that the fraction of 'fake' intensity in the observed 1 h wave in the Monte 
Carlo case is less than 5%, and thus negligible. 

In order to test the significance of the exotic wave, a second fit in mass bins was 
performed excluding the exotic wave from the wave set. A likelihood ratio test yields a 
log-likelihood difference of 47.3 between the two fits, averaged over a mass range of twice 
the experimental width around the resonance mass of the 7?! (1600). For a difference in 
the numbers of degrees of freedom of 4 this confirms the presence of the exotic wave in 
the wave set with a probability very close to unity. 

In conclusion, a partial wave analysis of COMPASS data from the diffractive dis- 
sociation of 190GeV/c 7T" on a Pb target into the 7r _ 7r~7r + final state at 0.1GeV 2 /c 2 < 
t' < 1.0GeV 2 /c 2 was performed. In addition to well-known qq' states, a significant natural 
parity exchange production of a spin-exotic wave with J PC = 1~ + decaying to pn is found, 
with an intensity of the resonant part corresponding to (1.7 ± 0.2)% of the total intensity 
in the mass-dependent fit. Its mass-dependent phase differences to the J PC = 2~ + and 
1 ++ waves are consistent with the highly debated 7r 1 (1600) meson. 

We gratefully acknowledge the support of the CERN management and staff as well 
as the skills and efforts of the technicians of the collaborating institutions. 
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